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Introduction 


Tropical rain forests contain a majority of the Earth's gene pool ( 
yet little is known of the spatial and temporal variation of these 
vegetation communities on a regional or synoptic basis. With the spread 
of man's activities the ability to monitor and forecast the response of 
these communities is becoming increasingly critical to preservation and 
management of such resources. 

This technical report represents the results cf research 
concerning two fundamental problems when considering the nonitoring of 
tropical vegetation succession. First, there is the 'shadowing' problem 
which is endemic to the use of Landsat in tropical areas. Secondly 
central to any monitoring system is the ability to model changes over 
space and through time for explanatory and predictive purposes. In the 
case of both problems the potentials of the field of spatial modeling 
were investigated. 

The Shadowing Problem and £ Spatial Model 

The Problem 

The problem of tropical natural resource inventory is fundamental 
to development of a monitoring capability. One of the major problems in 
the collection of data from Landsat for tropical areas is obstruction 
from clouds. The problem considered here is one where two statistically 
distinct classes occur with reference to an area where one homogeneous 
class of vegetation is known to exist. One class is homogeneous 
tropical rainforest and the other is known to be shadow. Variations in 
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data values which contribute to this dichotomy are posited as 
measurements subject influenced by variations at two scales, local and 
regional. There is a general pattern of spatial variation in both 
terrain and multispectral scanner (MSS) vaLues. In addition, there is a 
local pattern of variation unique to particular locations. It is 
therefore hypothesized that the shadowing effect is largely the result 
of small-scale spatial variations producing the responses measured by 
the values contained in the MSS data, thus causing the erroneous 
occurrence of two statistical classes where no observeable 
differentiation exits. The trend surface terrain model is used to 
decompose the spatial variation in elevation and MSS data. Since shadow 
is thought to be a function of the local component of spatial 
variation, the general component of spatial variation is used as input 
in the classification process. It is expected that after having 
filtered out the local spatial variation the two classes will become 
less distinct or not distinct. 

The Spatial Model 

The spatial model used here is generally known as the Trend 
Surface Analysis Model (TSA). It is used here as a means of 
quantitatively estimating the general spatial variation in the data 
set. In TSA the locations of the sample points provide information 
about the spatial variation of a variable. Each set of coordinates 
(U[i] and V[j]) describe the location of sample point k. The value of 
the observation (Z[k]) at point k is assumed to consist of 


2 



r 


contributions from two sources : 

REGIONAL or large-scale process which generates the broad pattern 
of spatial variation, and 

LOCAL or small-scale process operating locally. 

Fallowing Barringer et al (1980) and Robinson (1982), the model is 
specified as 

Z(U , V ) - f (U , V ) ♦ e 
i i i j ij 

where 

Z(U,V) ■ observed value at the location described by the 
coordinates U and V, 

f(U,V) - function defining the regional component, 

e ■ the error term which describes the local component 
since it is defined as Z(U,V) - f(U,V). 

The regional component function defined by fitting a power-series 
polynomial function of the general linear form 

n m i j 

Z(U,V) . E E b (U ,V ) 
k i«0 j»0 ij k 

which is estimated using the orthogonal polynomial routines presented 
in Mather( 1976) . In the case of Landsat data the measurement on Z are 
the data values on MSS bands, rfhile U and V are locations of their 
respective pixels. 

The Methodology 

The methodology was designed to evaluate the utility of a TSA 
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model in the reduction of misclassification of shadowed tropical rain 
forest. To accomplish this the following phases where were conducted in 
order of presentation : 


Phase 1 : Supervised classification of vegetation and shadow classes. 

Phase 2 : Perform linear discriminant analysis to test efficacy of 

supervised classification results and provide a benchmark. 


Phase 3 : Decomposition and separation of the regional and local 

spatial variations in spectral and terrain characteristics. 

Phase 4 : Using classifications of 1 and data generated in 3, a 

discriminant analysis was performed. 


Phase 5 : Comparison of the two discriminant analyses was conducted to 

evaluate the utility of the spatial modeling approach with 
regard to the supervised classification. 


Study Area 

The study area lies within a region of mountainous terrain located 
northeast of Mount Indie (approximately 12* 55' N , 121° 4' E) on 
Mindoro Island, Philippines. Primary cover is semi -deciduous rain 
forest comprised largely of dipterocarp spp. In addition, there are 
areas of (Ij_ cylindrica ) believed to be human induced. The topography 
is a dendritic drainage pattern. Thus, this is an ideal area for 
studying the shadowing problem. 


Cover Classification 

Supervised classification was conducted to identify the major land 
cover classes derived based on those derived from 1:50,000 vegetation 
overlays developed from aerial photography and ground data (Bruce, 
1979). Three major information classes were identified in a 128 by 128 
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pixel subscene. They were tropical rain forest, grassland, and shadow. 
For the conduct of the TSA modeling exercise, an area known to be 
tropical rain forest but shadowed and surrounded by rain forest. 

Sampling Design 

A random-systematic sample of elevation data and corresponding 
radiance count data was collected from within the study area. There 
were 89 pixels included in the sample. For each of the sampled pixels, 
elevation to the nearest 20 meter contour interval was interpreted from 
a 1:50,000 topographic map (AMS Series L501, sheet 3258 III, 1962). The 
elevation data was combined with the radiance count values taken from 
the four Landsat bands to form a data set which contains information on 
the relative location, elevation, and spectral characteristics of each 
sample pixel. 

Results of TSA Modeling 

Results of the TSA modeling of elevation and Landsat MSS data are 
presented below along with results of the discriminant analysis. 
Results of the discriminant analysis of the supervised classification 
are compared with results of the discriminant analysis using the TSA 
estimates as the data for cover class discrimination. 

The TSA Models 

A fifth-order polynomial trend surface was chosen as that order of 
surface which best describes the regional terrain variation within the 
sample area. The fifth-order model was chosen because it represents the 
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point of decreasing percent explanation added and the point where 
additional terms in the model provide less than modest increases in 
percent explanation (Table I). It should be noted that when one is 
concerned with describing a given pattern in terms of generalization at 
various scales, and relating those generalized patterns to the 
distribution of other variables, significance testing is unnecessary 
and inapplicable. 

A TSA model was estimated for each of the Landsat MSS bands. The 
order of surfaces chosen and the percent explained are reported in 
Table II. The percent explained refers to the percent variation in the 
original data accounted for by the regional component of the estimated 
surface. Thus, it is a measure of the relative importance of the 
regional component in explaining the spatial variations in the spectral 
responses. As can be seen irom Table II., the relative importance of 
the regional component in explaining variations in bands six and seven 
is small. The implication is that the local component in the spatial 
variation of spectral responses in these bands are primarily of a 
localized nature. 

The Discriminant Analyses 

Discriminant analysis of the supervised classification shows tha<: 
bands six and seven have the greatest discriminating power with all 
bands being present and significant in the linear discriminant function 
(LDF) (Table III, Univariate Analysis). The power of this discriminant 
function is exemplified by the classification results in which 83.15Z 
of the sample pixels were correctly classified. Since the initial 


6 



ORIGINAL PAGE 18 
OF POOR QUALITY 


classification was based on the MSS data, the significance and strength 
of this LDF is not unexpected. Thus, the shadow/f orest classification 
scheme is statistically significant when using the MSS radiance count 
data. 

When a discriminant analysis was performed with data generated by 
the TSA models, the results improved noticeably. That is the LDF is 
less able to separate shadowed and unshadowed forest. This is supported 
by an increased tendency toward mutual eauality (i.e., classification 
by random chance) of the four cells of the classification results table 
(Table IV). A fall in the eigenvalue and an increase in Wilks' Lambda 
are indicative of a decline in the discriminating power of this 
discriminant function. The all-groups stacked histogram (Figure 1) 
illustrates the variability and overlap present in this classification. 
Although the spectral terrain models for band 6 and 7 are those with 
teh weakest regional component in terma of explained variance, they 
maintain their dominant influence in the LDF (Table IV). It is 
important to keep in mind that their dominance in the LDF is 
conditional on the other variables being included in the estimation of 
the function. This is illustrated by the fact that when a stepwise 
discriminant analysis was performed, the only variables entered (in 
order of entrance were bands 7 and the regional terrain (elevation) 
component . 


7 


ORIGINAL PAGE IS 

OF POOR QUALITY 


Discussion 

This analysis supports cha view that shadowing effects may be due 
primarily to local variations in the spectral responses. The 
significance of this result is that they can be compensated for through 
the decomposition of the spatial variation in both elevation and MSS 
data. Use of the trend surface model to estimate both elevation and 
spectral terrain surface aa a. posteriori inputs in the classification 
process leads to an improvement in classification accuracy for 
vegetation cover of this type. The improvement might be increased 
further through use of weighted linear components of the MSS and 
elevation terrain data or the specification of a power-series ocher 
than the straightforward polynomial expansion. 

Results of this portion of the research alsj suggest that the 
spatial patterns depicted by the “IS data reflect the measurement of 
responses to spatial processes sewing at several scales. Thus, 
continued research on the development and application of spatial 
filtering models in the analysis of Landsat data is suggested. Perhaps 
a more broadly significant implication is the use of dynamic spatial 
models will be critical in the study of tropical rain forest succession 
using Landsat. Therefore, the latter portion of this research 
investigates the field of dynamic spatial cover class models as they 
relat? to succession dynamics. 
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Vegetation Succession, Land sat , and Dynamic Spatial Models 

An Investigation 

Here are presented important issues which must be carefully and 
thoughtfully considered if remotely sensed data are to be properly used 
in the analysis of modeling of spatial succession processes, 
particulary in the tropical rain forest context. Alinough raised in 
the context of tropical vegetatiou cover dynamics, many of these issues 
have broad applicability to other research using Landsat data to 
monitor changing vegetation patterns (e.g. see Robinson et al., 1982). 

Remote Sensing of Vegetation 

Efforts in the field of satellite remote sensing of vegetation 
have traditionally centered on questions of signature identification 
and extension, cover inventory accuracy, and change detection and 
monitoring (e.g., Anderson, 1977 ; Bannert, 1980 ; Omakupt and 

Vunpiyarat, l! ' ; beeves, 1975 ). Although the remote sensing 

literature contains many references to the spatial and temporal aspects 
of the science (e.g., Lillesand and Kiefer, 1979 ; Swain and Davis, 
1978). Approaches incorporating the spatial and temporal dimensions are 
typically not of a formal nature. Formal models of cover dynamics are 
often specified in aspatial form even though satellite data are clearly 
3patial in nature. Also, design of satellite spatial sampling 
instrumentation is largely a function of scanner technology and image 
resolution considerations. Little consideration seems to have 
been given to the unique statistical properties of spatial samples 


13 



ORIGINAL PAGE IS 
OF POOR QUALITY 


(e.g., Baccy and Sammons, 1978 ; Bellhouse, 1977 ; Curtis and 
McIntosh, 1950). 

Models of Succession 

The theory of ecological succession has been developed in Che 
papers of Clements (1916, 1936), Gleason (1926), Whittaker (1953), Odum 
(1969), and Drury and Nisbet (1973). Although the term was coined in 
the early 1800's, Clements is generally regarded as the first to offer 
a theory of the succession process. He argued that "The developmental 
study of vegetation necessarily rests upon the assumption that the unit 
or climax formation is an organic entity... As an organism the 
formation arises, grows, matures and dies" (Clements, 1916). In 
contrast to Clements (1916), Gleason (1926) took the view that the 
successional processes of a plant community grow out of the processes 
which its individual members undergo: "The sole conclusion we can 
draw... is that the vegetation of an area is... the resultant of t wo 
factors, the fluctuating and fortuitous immigration of plants and an 
equally fluctuating and variable environment" (Gleason, 1926). 

Another view of succession is that presented by 
Whittaker (1953). He believed succession to be due to interactions 
among a series of environmental gradients: "Vegetation is conceived 
as,... a pattern of populations, variously related to one another, 
corresponding to the pattern of environmental gradients" (Whittaker, 
1953, pp. 59-60). 
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Miles (1979) describes Odum's (1969) model as neo-Clementsian 
because of its emphasis on properties emergent at the community level. 
Odum (1969) saw the "strategy" of succession as being action by the 
community to bring about an increased control of, or homeostasis with, 
the physical environment to achieve maximum protection from its 
perturbations. Drury and Nisbet (1973) argued that succession grows 
out of differential colonizing ability, growth and survival of species 
adapted to growth under different environments. 

Some general observations may be made on these models. 
First, they are all informal models (van Hulst, 1978). Second, while 
all of them imply the existence of temporal and spatial features in the 
process, none of them adopts a formal notion of either. Neither do 
they consider the essential differences which exist between spatial and 
temporal processes. Indeed, Drury and Nisbet (1973) note that most 
descriptions of long-term changes in vegetation are based on the 
observation of spatial sequences. Third, the models tend to shift 
increasingly from a deterministic one toward a stochastic one (see van 
Hulst, 1974). Finally, the neo-Clementsian arguments have been 
essentially demolished by Drury and Nisbet so that some trend from a 
holistic toward an atomistic view of succession is evident. 

Beeftink (1979) presents several reasons for increased 
interest in vegetation dynamics on the part of European ecologists. 
First, interest in descriptive research is falling off for a variety of 
reasons. Next, vegetation conservation efforts demand a better 
understanding of vegetational and environmer. *al dynamics. Last, 
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advances in numerical methods and instrumental equipment present the 
opportunity of developing dynamic models utilizing large data sets. 
He recognizes several points of departure for future research in 
successional dynamics. One such question is that of the degree to 
which detection of similarities between space-time patterns in 
biological and in environmental variables provide insight into causal 
relationships among these variables. The connection between diversity 
and stability concepts, and space-time patterns is central to such 
efforts. 

Allen (1975) observes that in the modeling of population 
stability over time, the usual approach ignores population variation 
across space. He comments on the strangeness of this and warns that 
time-dependent models can be radically affected by extending them into 
the spatial domain. Therefore, the theoretical and experimental 
investment of such models is certainly justified. Moreover, we might 
seriously question whether stability statements based on the original 
models can be trusted. Yarranton and Morrison (1974) investigate the 
variation in spatial pattern over time in primary successional 
sequences. They present a nucleation model which they argue 
establishes a relationship between successional mechanism and 
variations in spatial pattern. Hilbourne (1979) points to the 
increasing realization that spatial phenomena are an important aspect 
influencing the dynamics of many ecosystems. He argues that an 
understanding of the formation of spatial heterogenieties is crucial in 
understanding the temporal dynamics of ecosystems. Finally, Miles 
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(1979) warns against the interpretation of existing spatial variation 
in vegetation as representing different phases in a temporal sequence. 

It is apparent that there is considerable interest among 
ecologists in successional processes as they occur in both time and 
space. Similar interest is in evidence among those of the remote 
sensing field and biogeographers (e.g., Miller, et al. t 1978). If 
these interests are to be addressed, two important issues must be 
addressed. The first complex of issues concern those the properties of 
static spatial models. Secondly there are the implications of spatial 
sampling methods for the study of succeessional processes. Although 
such static models permit little insight into the dynamics of 
succession, they form an important class of descriptive models and In 
their descriptive role, they may aid in hypothesis formation for 
dynamic models. 

Static Spatial Models 

Hera are presented those spatial models of a static nature 
which are important in developing models of tropical rain forest 
succession using Landsat data. It is common to regard the distribution 
of a set of points in geographical space at some given time as the 
outcome of a point process. While the actual generation of such a 
distribution is presumed to be dynamic, the resulting distribution is 
not. More importantly, nor are the underlying spatial models. 
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Bartlett (1975) organizes point process models into two broad 
classes: 

A. l Clustering models 

.2 Doubly stochastic models 

B. Contagion or inhibitory models 

All are based on the Poisson distribution. More complex point 
patterns may be described by mixed Poisson models which represent a 
mixing, in various ways, of Poisson process rules with those of other 
statistical distributions (e.g., Poisson-Poisson, Poisson-binomial). 

Data for such models are usually obtained by the quadrat method. 
First, a planar region is exhaustively partitioned into a series of 
non-overlapping, generally square cells. Next, the numbers of points 
falling in each cell are determined. Finally, the degree of randomness 
of the distribution is examined through the fitting of the 
distributions described above which describe the degree of clustering 
or regularity of the point set. 

Rogers (1974) and Curtis and McIntosh (1950) point out some 
of the difficulties associated with the quadrat method. Results depend 
upon the shape, size and placement of quadrats as well as on the 
density of the points being sampled. They propose rules for the 
determination of optimal quadrat size, as does Grieg-Smith (1964). The 
relative effects of applying alternative optimal size rules to quadrat 
size selection should be of interest to those engaged in research 
involving spatial data collection. 
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Rogers (1974) also discusses inferential problems which must 
be considered in testing the goodness-of-fit of point process 
statistics to non-random (i.e., either highly clustered or highly 
regular) spatial data. The chi-square test, a statistic commonly used 
in such applications, is noted to have a number of shortcomings in this 
respect. Among these are: lack of sensitivity, influence of data 
grouping schemes on the outcome of the test, minimal class size 
constraints on test effectiveness in the face of a lack of consensus on 
a correct value for minimal size and loss of test power connected with 
grouping to satisfy minimal size constraints. Rogers also notes that 
similar difficulties attend the fitting of bivariate distributions. A 
far more serious problem of inference in that statistical agreement 
with a mathematical model is far from sufficient evidence that it did 
so arise. 

The statistics of nearest neighbor models also raise issues 
for those using them for spatial inference. It has been noted (Rogers, 
1974) that density dependence, derivation of probability distributions 
for non-random point patterns and the effects of 'zero-distance' 
nearest neighbor configurations are all questions of interest in this 
area. 

Description of spatial arrays of points is closely tied to 
the effect of process scale on scale of sample and the correct design 
of spatial sampling schemes. Clearly, the relative spatial extents of 
a spatial process and of a spatial sampling scheme designed to provide 
data for inference on that process will have an effect on the outcome 
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of such inference. Also, Kooijman (1980) states that the ease of 
inference on point processes decreases with decreasing spatial scale. 
Finally, the optimal shape is considered to be a rectangle (Or loci, 
1978). 

Haggett, et al. (1977) examines the relative efficiencies of 
various spatial sampling schemes. In that review, Berry's (1962) 
determination of the high relative efficiency of stratified systematic 
unaligned sampling schemes over either stratified random samples or 
simple random samples is cited. Further, the superiority of stratified 
sampling plans in situations where spatial autocorrelation may be 
present is noted. Haggett (1963) also found transect sampling methods 
to be more accurate than either point or quadrat sampling in measuring 
the extent of cover by vegetation in an area of England. 

Rogers (1974) looked at stratified random sampling in a 
spatial context and Bellhouse (1977) derives a minimal average variance 
criterion for sampling from a two-dimensional finite population. He 
also concludes that no global optimal design exists although locally 
optimal designs do exist in terms of certain subclasses of spatial 
sampling scheme. A problem, until quite recently, in spatial sampling 
schemes was that of the asymptotic distribution of the sample. This 
issue was resolved with Smith's (1980) publication of a Spatial Central 
Limit Theorem. Although the theorem gives little guidance on the size 
of the sample needed to realize asymptotic normality, it clearly shows 
the existence of such a property for a spatial sample of sufficient 
size and meeting certain other conditions. 
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Issues In the Analysis of Spatial Data 

Although Landsat data are quite obviously spatial data, for the 
purposes of classification and modeling, reliance has been placed 
primarily upon classical statistical models (e.g., Jensen, 1978; 
Strahler, et al., 1980). Tubbs and Coberly (1978) pointed out that most 
investigators have tended to assume that pixels are independent 
observations. The fact that neighboring pixels are not separate 
independent observations is at the core of the concept of spatial 
modeling. 

There are several basic properties of spatial data which make 
their analysis us^ng classical statistical models difficult. Spatial 
statlonarity is an implicit assumption in all applications of Landsat 
data using classical statistical methods. In general, spatial data do 
not possess the property of spatial statlonarity. Granger (1969) 
describes statlonarity as an assumption that the relationship between 
values of the processes generating the data is the same for every pair 
of points whose relative positions are the same. Cliff and Ord (1975) 
have presented a more formal definition of spatial statlonarity. First, 
it is assumed that the data collected correspond to a finite set of 
locations. In Landsat data the locations are analogous to the pixels in 
a Landsat scene or subscene. Also, in the Landsat case, the 'spatial 
aggregates' are areas not points. Next, let j denote the set of 
locations and j ■ l,2,...,n where n is the total number of locations. 
There is the observed variate value, y or Y . There is information 

j J 
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concerning location j , e.g. row and column location of a pixel. 
Suppose that Y can be decomposed into a stochastic component X and 

J J 

a deterministic component m such that: 

J 

E[ X ] - 0 or E[ Y ] - n> . 

J J J 

Definition 1. X is said to describe a spatially stationary 
j 

process in the wide sense, to be weakly spatially stationery, if: 

E( X ,X ] - or (j,j’) (1) 

i J' 

depend, only upon the relative position of locations j and j'. 

Definition 2. Suppose X satisfies equation (1). In addition, 

J 

if the correlation between X and X depends only upon the distance 

J i' 

between their locations, and not upon the orientation between j and 

J' , then the process is weakly isotropic. 

Definition 3. If the joint distribution of the X depends onlv 

J 

upon the relative positions of the locations, then the spatial process 
is strictly stationary. 

Definition 4. The process is strictly isotropic if it is both 
strictly stationary and direction invariant see Definition 2). 

In areas such as time-series analysis definitions 1 and 3 will 
suffice, but in spatial analysis definitions 2 and 4 are often required 
(Cliff and Ord, 1975; Haggett, et al., 1977). Generally, the assumption 
of spatial stationarity is difficult to sustain. It is worth noting 
that the concepts of spatial stationarity and spatial dependence are 
related. Haggett, et al., (1977) point out that if the marginal 
distributions of X are identical, then non-stationarity implies that 

j 
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the data ara spatially dapandant. 

Spatial dapandanca among elements of spatial data sats is tha 

rule rathar than tha exception (Toblar, 1970; Gould, 1970). To 

furthar illuatrata this, supposa that data hava been collactad for a 

sat of j ■ l,2,3,...,n araaa (or pixals) and that tha variata for tha 

j-th location is X . Statistical aathods ara traditionally davalopad 

J 

from assuaptlons which usually stata something Ilka - 
” let X ( j-1,2,. . .,n) be n independent, identically distributed 

J 

variates” (Cliff and Ord, 1975a in Haggett, at al., 1977). 

Spatially located data generally exhibit systematic spatial 
variation or spatial autocorrelation (Toblar, 1970; Gould, 1970; Cliff 
and Ord, 1975). Tha affect of spatially autocorrelated data on tests 
of inference ara in soma cases wall known. Standard applications of t 
and F statistics for comparison of means or construction of confidence 
intervals require spatial independence. Tha same assumption is 
necessary in regression analysis if the OLS (Ordinary Least Squares) 
estimators are to be BLU (Best Linear Unbiased). 

Many common applications of a variety of statistical models in 
remote sensing use as the variance estimate: 

- I (y - 7 )* / (n - 1) 

J-l J 

Cliff and Ord (1975) note that when positive spatial 
autocorrelation (e.g., a clustering process) is present, the estimator 
9 * will be biased downwards. This bias leads often times to an 
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overstateaent of the significance of the results. This Is an lssus not 
widely recognized in the applicecion of statistical aodels co Landsat 
data. 

It has also been shown (Brandsaa and Ketellapper, 1979) that che 
choice of a procedure for eatiaating test statistics for spatial 
autocorrelation oust be aade carefully since the power, robustness and 
obtainable significance levels say vary considerably aaong estiaatlon 
procedures. 

One of the aore coaaon statistical aodels used in the analysis of 
Landsat data for cover class studies is principal coaponents analysis 
(PCA). Lebert's (1969) approach is briefly presented below as an 
ezaaple of how spatial dependence aay be incorporated into this coaaon 
aethod. Suppose data are collected on p variates, say Landsat MSS 
bands 4,5,6, & 7 , for each of n areas ( pixels ) . This yields a (n x p) 
data aatrlx X. In addition, a aatrix M is constructed where, 

a • 1 if the 1-th and J-th pixels are contiguous, 

ij 

■ 0 otherwise. 

M ■ M**a defines the nuaber of paths , a links in length, between 
a 

each pair of pixels. This definition Includes redundant paths. The 
saaple a-lag covariance for the k and 1 variates (e.g. bands 4 and 
7) is in aatrlx fora: 

T 

C - (l/n ) X (N - M ) X (a - 0.1,...) 

a a a a 
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2 

where n is the number of paths of length a , and N ■ diag(M ) It may 
a a a 

aid in understanding to nota that C la the usual sample covariance 

0 

matrix. Lsbart (1969) argues chat the use of C other than C will 

a 0 

reduce the effects of spatial dependence upon the analysis. 

The example dealing with PCA begs the question of whether one 
should: (1) Incorporate spatial dependence into the model explicitly, 

or (2) remove it before commencing the main analysis. The choice of 
approaches depends upon the emphasis required by the nature of the 
problem to be solved. However, from the perspective of spatial 
modeling one aay argue that corrections used to remove spatial 
dependence represent 

N . . . a throwing out of the uaby and keeping the bathwater " 

(Gould, 1970). 

Another fundamental problem of spatial data has to do with aggregation 
problems. The results of moat any analysis depend on the manner in 
which the geographical study area la partitioned for data collection 
and processing purposes. We can use the example of Yul« and Kendall 
(1957) in whici the value of their computed correlation coefficient (r) 
between yields per acre of wtieat and potatoes varied rather widely 
depending upon the 'spatial aggregates' used. This problem is directly 
comparable with the same methodology reported by Robinson et al (1982) 
in a study of desertification modeling using Landsat data. In contrast 
to Yule and Kendall (1957), Robinson et al (1982) does not appear to 
have been aware of the role of spatial aggregation problems when 
partitioning the study area. Even so, his study does represent a 
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methodology which recognizes the spatial variation in variate 
relationships. In other words, Robinson et al(1982) study recognize in 
a rudimentary way the fact of spatial non-stationarity. 

Bearing in mind the dual issues of spatial stationarity and 
autocorrelation, the analysis of spatial data usually has as its 
objective either: (1) data smoothing, (2) data interpolation, or (3) 
modeling. These objectives are not necessarily mutually exclusive. For 
example, the preprocessing procedures used to destripe Landsat data 
introduce, in an obvious manner, strong components of spatial 
autocorrelation before any modeling has been done. This is also the 
case in data interpolation for bad or missing Landsat data. 


Spatial Modeling 

The natural framework for the modeling of cover class dynamics is 
the spatial-temporal process diice we need to capture both components 
through temporal and spatial dependencies. However, before spatio- 
temporal models can be developed, spatial models must be considered. 
Purely spatial processes are philosophically highly problematic. Their 
study is of interest for three reasons (Bennett and Chorley, 1978): 


(1) They may, under special conditions, actually exist. Spatial 
processes actually exist when the system under study has reached 
equilibrium. 

(2) They may seem to exist to the observer or analyst. This usually 
occurs when the space-time dat.i are available with very wide time 
separation which does not allow recognition of reaction, relaxation, 
and lagged dynamics. 

(3) There are a number of important spatial problems such as spatial 
extrapolation, interpolation, and selective problems of statistical 
inference . 
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It is mathematically straightforward to develop representations 
of purely spatial models. Two spatial processes which have received a 
great deal of attention are the unilateral and multilateral models 
defined on the basis of spatial lag as: 

unilateral 

Y - a (Y + Y )+ e 

tij ti(j-k) t(i-l)j tij 

multilateral 

Y ■ a (Y + Y ♦ Y + Y )+ e 

tij t(i-k) j t(i+k)j ti(j-l) ti(j+l) tij 

k,l ■ 1,2,3, ... 

where a Cartesian lattice is assumed. The parameter a provides an 
indication of spatial stationarity where in each case the process is 
stationary if for the 

unilateral case: 

lal < 0.5 

multilateral case: 

lal < 0.25 

(Haining, 1978). 

Multilateral models have been used to describe interplant 
competition in ecology (Bartlett, 1975; Mead, 1971). The important 
difference between multilateral and unilateral models is that 
multilateral models are space-like whereas unilateral models perform in 
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a way analogous Co an ordered time-series. In ocher words, in a 
unilateral model there is a distinctive ordering of events as they 
occur across the plane. An event initiated at one location spreads out 
in only one direction (Bennett, 1980). 

There have been few formal models of Landsat data which have 
recognized the problem of spatial dependence. Notable attempts are 
those of Tubbs and Coberly (1978), Craig (1979), and Craig and Labovitz 
(1980). Without exception, attempts to date have depended heavily upon 
some version of the unilateral model. Hence, it is only natural that 
they should rely totally upon the Box-Jenkins (1970) formulations of 
the AR (autoregressive) and ARIMA (autoregressive integrated moving 
average) models. In contrast to these formal modeling attempts stand 
the less formal interpolation or filtering functions commonly used in 
the preprocessing of Landsat data. These filtering functions are, in 
the main, based on the same concept as the multilateral model. In 
fact, the destriping algorithm is essentially a version of a spatial 
moving average of a multilateral nature (Taranik, 1978). This implies 
that those using such preprocessed data should take multilateral 
spatial autocorrelation into account for the data they are working with 
are very much related to the data values in surrounding pixels. 

The approaches taken by Craig & Labovitz (1980), Craig (1979), 
and Tubbs 4 Coberly (1978) underscore the relationship between the 
spatial autoregressive and the temporal (unilateral) autoregressive 
schemes. Let us consider the time-series autoregressive model of first- 
order where X(t) is the variate measured at time t (t ■ ... , -1,0,1, 
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... ). Three equivalent definitions follow which maintain the 

condition of stationarity by constraining !p! < 1.0 (Cliff, 1980 ). 

The simultaneous autoregressive model is: 

X ■ pX + e where, 

t t-1 t 

E(X ) » 0 ; E(e )-0 ; var(e ) ■ o’ ; cov(e , e ) » 0 ; 
t t t t t-1 

cov(X , e ) * 0 for s > 0 . 
t-s t 

The conditional autoregressive model is: 

EC X IX , X ,...)- pX 
t t-1 t-2 t-1 

and 

Var(X IX , X ,...)- <r* 
t t-1 t-2 

The covariance model is: 

E(X ) - 0 and cov(X , X ) 
t t t-s 

where, er* ■ c r 1 / ( 1 - p* ) 

x 

(Cliff, 1980). 

Whittle (1954) and Granger (1969) have shown that when a first- 
order autoregressive model such as: 

X - pX + e 
t t-1 t 


1st 


» cr p 

X 
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is expanded co Che two-dimensional , multilateral form: 

X • aX + bX + cX + dX +e , 
ij i-l.J i+l.J i.j-1 i.j+1 ij 

considerable problems are encountered. The two major problems are (1) 

parameter estimation becomes computationally difficult and (2) on the 

basis of the same set of spatial autocorrelations, it is impossible to 

discriminate between a number of alternative models (Bassett, 1972). 

This problem of model identification or specification will be treated 

in greater detail in the discussion of dynamic spatial models. 

The autoregressive and moving average schemes assume spatial 
data in discrete space. However, digital Landsat data are obtained 
through an analog-to-digital conversion process. In this process the 
[continuous] signal is sampled frequently enough so that the digital 
representation of the signal will reproduce the information content of 
the signal (Swain and Davis, 1978). This characteristic of Landsat 
data suggests the possibility that models based on the discrete-space 
assumption may be inappropriate and those based on continuous, but 
sampled, spatial data are more justified. If this is the case, Cliff 
(1980) suggests that it is better to specify the covariance structure 
directly rather that attempt to formulate a linear dependence model. 
The reasoning is that if a linear scheme is specified at one spatial 
scale, there would need to be a scheme which is non-linear in p at 
any different spatial scale. 


30 



ORIGINAL PAGE IS 
OF POOR QUALITY 


We may consider what the general formulation of such a covariance 

structure might be. Consider a process X(u) defined at: 

location u ■ ( u , u ) 

1 2 


with 


mean E [ X(u) ] ■ u 

and 


covariance structure 

cov[ X(u),X(v) ] • or* c (u , v) 
where c(u,u) ■ 1 for all u . 

Various forms of c(u,v) have been suggested usually depending on 
a distance formulation such as: 

d a [( u - v)* + (u - v)* ]$ , 

rewriting c(u t v) as c(d) and keeping c(0) a 1.0 

This is only one of several alternative approaches to 
specification of the form of the covariance structure (see Cliff, 
1980 ). 


Temporal Series and Cover Class Dynam-i cs 

When speaking of cover class dynamics one is generally interested 

in process. Process is essentially change over time. Time-based systems 

possess convenient properties not encountered in spatial systems . The 

most important properties are asymmetry and transitivity. 

Mathematically these properties can be defined for events E(i) and time 

operators T by 
i 
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Asymmetry : 

given E - T (E ) , then E / T (E ) 

112 2 11 

Transitivity : 

given E » T (E ) and E - T (E ), 

1 1 2 2 2 3 

then E « T (E ) . 

1 3 3 

In other words, the whole set of events is uniquely connected or 

ordered. One of the simpler models suggested for use in the study of 

land cover change using remote sensing data (Bell, 1974; Strahler, et 

al., 1980) and used extensively in the study of plant succession 

(Slatyer, 1977) is the Markov Chain Model. The Markov Model is 

concerned primarily with changes of 'state' over time. The transition # 
matrix forms the basis for the definition of the model where the 
transition probability matrix P defines the probability of cover 
class i changing to to cover class j from time t to time t+l . A 

transition probability is the probability that a given spatial unit S 

in cover class i at time t will be in cover class j at time t+l. In 

notation form, elements of the transition probability matrix may be 
defined as : 


t,t+l 

P - Prob [ S - j ! S - j ] 
ij t+l t 

(Robinson, 1980). 

Bell (1974) suggested its use in the study of land use change 
data derived from aerial photographs and Strahler, et al. (1980) have 
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used it as a means of analyzing change in forest cover classes. From 
the perspective of studying cover class dynamics as a process of plant 
succession, Markov Models have the following advantages : (1) they are 
easy to derive from successional data; (2) they do not require deep 
insight into the mechanisms of dynamic changes, but it can help to 
pinpoint areas where such insight would be valuable; (3) the transition 
matrix summarizes the essential parameters of dynamic change in a 
system in a succinct, concise manner; (4) the computational 
requirements are modest. On the other hand, there are several serioua 
disadvantages such as : 

(1) the lack of dependence on functional mechanisms may reduce 
its .appeal to those less empirically-oriented scientists (Slatyer, 
1977); 

(2) departure from assumptions of stationarity, or from 
assumptions of constant, linear transfer rates (Bell and Hinojosa, 
1977); 

(3) validation of the model depends upon predictions of system 
behavior over relatively long periods of time and is therefore 
difficult in many cases (Robinson, 1980). 

There are two additional problems of the Markov Model which are 
associated with the use of Landsat data for modeling cover class 
dynamics. These are: classification and spatial dependence. In order 
to use the Markov Model to study cover class dynamics, some 
classification scheme is needed which separates successional 
communities into definable categories such as cover-states (Slatyer, 
1977). The problem of classification remains an elusive problem in the 
study of land cover change (Stow, et al., 1980). The second problem 
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has received less ac tendon in the remote sensing literature. The 
problem of spatial dependence was recognized in the study by Bell and 
Hinojosa (1977) as being that the Markov model is intrinsically 
aspatial. Thus, there remains the problem of finding a variation [of 
the Markov model] that can take account of this dependence explicitly. 

Uncritical use of the Markov model results in the same problems 
as in the uncritical use of other, more general, models of time-like 
systems. Namely, the problems of stationarity and autocorrelation 
affect both the Markov model and the general linear stochastic models. 
Autocorrelation (in the non-spatlal sense) refers to dependence (i.e., 
correlation) between the observed value of a variable at some time t 
and its observed value at some later time t+1. As such, it represents 

a violation of one of the assumptions of the Gaussian model. 
Stationarity refers to constancy in the moments of a time series. A 
series is referred to as being first-order stationary if its mean is 
constant from one time period to the next. Second-order stationarity 
is similarly defined with respect to variance. If a series is second- 
order stationary, the implication is that the covariance of two terms 
of the series is dependent only upon time difference between the points 
and not upon the particular points being examined (Cliff, et al., 
1975). That is: 

cov[X(t),X(t-k)] - E[(X(t) - u )(X(t - k) - u)]. 

The above reference quotes Granger (1969) to the effect that 'The 
assumption of stationarity essentially says that the law that generates 
the data is constant over time'. 
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An Introduction to Dynamic Spatial Models 

Despite a growing interest in space-time modeling, 
especially in the remote sensing community, there is a relatively scant 
body of literature on the topic. This section presents three general 
spatial-temporal models, they will be the STAR (Space-Time 
Autoregressive Model), STMA (Space-Time Moving Average Model), and the 
STARMA (Space-Time Autoregressive-Moving Average Model) models. 

As is evident from the types of models which shall be discussed, 
this section is based on the modeling of Landsat data for 
classification, not on modeling classified Landsat data following an 
aspatial classification exercise. This approach is suggested by the 
work of Tubbs and Coberly (1978), Craig (1979), and Craig and Labovitz 
(1980). Furthermore, concern with modeling Landsat data prior to the 
classification exercise has offered promising results in the modeling 
of tropical rain forest succession (Barringer, et al., 1980), and in 
the monitoring of desertification (Robinson et al, 1982). In addition, 
this approach may been seen to act to 'filter' out certain kinds of 
confounding influences which arise from sources noted in detail by 
Craig and Labovitz (1980), and Barringer and Robinson (1981). 

To begin the discussion of these models some common definitions 
are presented. First is the definition of the "spatial lag operator". 
For a spatial lag operator L and data value y(it) that is associated 
with pixel i for time t, the relationship is: 
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and 


0 

L y - y 

It it 



9 

Ly ■ E w y ,s>0, 
it J ij J.t 

where 9 denotes the spatial lag which 19 9 step9( pixels) away from 
pixel i ( and the summation is over the j pixels at spatial lag s from 
1. Using this common notation all the above models are briefly 
presented . 

STAR Model 

Tobler (1967, 1970) used STAR models to estimate the linear 

spatial transfer functions (linear operator) which best transforms a 
map at time t into that at t+1. This is, of course, a common problem in 
the use of Land sat data to model cover class dynamics. The STAR model 
is of the fora: 

1 m s 

y - E E a L y + e . 

it s»0 k«l sk i,t-k it 


where k denotes temporal lag k from t, while 1 and m denote the maximum 
number (order) of the spatial -temporal lags over which the summations 
are conducted (Haggett, et al., 1978). 
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STMA Model 

This a space-time moving average model of the general form: 

1 p s 

y - - E EcLe 4-e. 

it a -0 k«l sk i,t-k it 

STARMA Model 

Combining the STAR and STMA models we obtain the space-time 

autoregressive and ooving average model where: 

lms 1 p s 

y ■ K E a L y - E EcLe +e 

it s-0 k-1 sk i,t-k s-0 k-0 sk i,t-k it 

Model Specification and Identification 

In the stochastic modeling of successional processes, ecological 
theory forms the basis for model specification. Box and Jenkins, 
(1976) point out that where theory is inadequate or non-existent, a 
class of models may be specified. However, a given model may imply 
more than one possible structure. 

Identification forms the second stage in the process of 
statistical inference and is the logical process whereby the model 
specified in the first stage is related to a basic structure. Correct 
model specification, based on ecological theory, implies structural 
identification in terms of the model. Moreover, the identification 
process works to reduce any class of models which may have been 
specified initially to a more parsimonious sub-class. Haggett, et al., 
(1977) note that two questions may be posed which together summarize 
the model specification problem. These two questions are: 
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(1) Which nodal should be selected from che basic STAR, STMA, and 
STARMA formulations prtaantad abova ? 

(2) Having choaan a particular aodal, what ordara of tenporal and 
spatial lags should ba includad in tha nodal ? 

If acologic thaory is strong anough to provlda, a priori , answars 
to thaaa two quastiona, than tha problaa of nodal apacification is 
aasily solvad. Howavar, more often than not, these bodies of thaory 
offer little guidance in naking such fine distinctions anong modeling 
forna. 

Tha identification of an appropriate aodal requires not only 
coaputation of tha tine-space correlogran (i.e., a display of tha 
serial correlation structure anong variate values for various spatial 
and taaporal lags, but also of tl time-space partial correlogran which 
displays partial correlation relationships between and anong the 
Independent variables. Tha partlals indicate the degree of correlation 
between the data values in the pixels at tine t and the values in those 
pixels which are k tonporal and s spatial lags away with the effect of 
all other spatially and tanporally lagged variables held constant. 
Each of the nodels presented above has been shown to have certain 
general characteristic space-tine correlations and partial 
correlations. They are: 

(1) For the STAR (l,a) process of order 1 in space and a in tine, 
the autocorrelations should decay approxinataly exponentially in tine 
and space. The partial correlations should become approximately zero 
after lag 1 in space and m in time. 


38 



ORIGINAL PAGE 18 

OF POOR QUALITY 


(2) For the STMA (g,p) process of order g in apace and p in time, 
the autocorrelations approach zero afcar lag g In apaca and p in time. 
Partial corralationa decay approximately exponentially in tlaa and 

apaca. 

(3) For cha STARMA (l,m,g,p) procaaaea diaplay corralationa and 
partial corralationa which tail off aa a aixtura of axponantial curvaa 
and daapad aina wavaa. Tha corralationa axhibit thla tailing-off aftar 
tha firat g-1 spatial laga and p-a laga in tima Tha partial 
corralationa axhibit auch bahavlor aftar tha firat 1-g laga in apaca 
and a-p laga in time (Martin and Oeppen, 1975). 

In ordar corractly to spacify and ldantify a apatial nodal of 
succaaaion froa ona of tha claasas outlinad abova, thought auat ba 
given to thraa aata of relationships: 

(1) For any given region in space and for any two points in 
time, what effect will tha apaclaa collection present on 
the site at tiaa t have on tha apaclaa collection present 
at tiaa t + k ? (Temporal autocorrelation affect). 

(2) For any given tiaa and for any two regions in apace, what 
affect will the apaclaa collection present at site i have on 
tha species collection at site j — and vlca versa? (Spatial 
autocorrelation affect). 

(3) For any regions i at time t and j at tiaa t ♦ k, what effect 
will tha species collection on site i have on that at site 

j — and similarly for tha effects of site j at tiaa t on 
1 at tiaa t + k ? (Space-time autocorrelation). 

In tha study of cover class dynaalcs aa a successlonal process 
such models may prove lass than illuminating if ona does not realize 
that they era based on tha assumption that the structural relationships 
between spatial aggregates may vary over space and time — in 
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particular, patterns of decline and cutoff for the various 
autocorrelation effects will be of considerable importance in the model 
specification. For example, it may be unrealistic to assume that the 
form of the relationship between two cover classes at their mutual 
boundary is the same as it i« in an area of homogeneous cover located 
away from that boundary. Haining (1978) encountered the related 
problem of bias associated with the estimation of certain stochastic 
lattice models. It was observed that such edge effects are dependent 
both upon the ratio of border to non-border cells (e.g., pixels) and 
upon the values of border cells relative to the expected values of non- 
border cells. Lattice regularity was also found to have a strong effect 
on process outcome. It is to be expected that effects such as those 
alluded to here would be particularly troublesome when studying such 
spatial systems over time. One sophisticated approach to this problem 
is the construction of varying parameter space-time models (Haggett, et 
al., 1977). 

Concluding Summary 

There are two major outcomes of this research. First, the use of 
the TSA model shows that the spatial modeling of radiance values can 
provide a useful approach to one the problems in monitoring tropical 
rain forest succession. It also raised the point that shadowing effects 
may be due primarily to local variations in the spectral responses. 
The significance of this result is that they can be compensated for 
through the decomposition of the spatial variation in both elevation 
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and MSS data. Use of the trend surface model to estimate both elevation 
and spectral terrain surface as £ posteriori inputs in the 
classification process leads to an improvement in classification 
accuracy for vegetation cover of this type. Perhaps more significant 
are the results suggesting that spatial patterns depicted by the MSS 
data reflect the measurement of responses to spatial processes acting 
at several scales. The results of the TSA modeling phase suggested the 
investigation of spatial models as useful in the development of a 
capability in monitoring tropical vegetation changes. 

The most broadly significant result of this research is the 
identification of dynamic spatial models critical to the study of 
tropical rain forest succession using Landsat. A major theme has been 
the spatiallity of Landsat data, of the process which it represents and 
of the models which may be used to emulate successional processes. 
Thus, its effect on data collection, model construction and statistical 
inference were identified. An effort was made to suggest the promise 
of dealing formally with space, as well as some of the pitfalls which 
may await the unwary. It is strongly felt that the issues raised here 
must be clearly considered and resolved if future models based on 
Landsat data are accurately to reflect the space-time processes they 
are meant to emulate. 


41 


Bibliography 


Allen, J.C., 1975, "Mathematical Models of Species Interactions 
in Time and Space", Amer. Wat. , 109(967): 319 - 342. 

Anderson, J.R., 1977, "Land Use and Land Cover Change: A 
Framework for Monitoring", Jrnl. Res. USGS , 5(2): 143 - 153. 

Bannert, D., 1980, "An Operational System to Monitor 
Desertification Using Landsat Data", Proc. 14th Inti. Symp. Rem. Sens. 
Envt. . Envt'l. Res. Inst. Mich. Ann Arbor. 

Barringer, T. and V.B. Robinson. 1981. "Stochastic Models of Cover 
Class Dynamics" . Proceedings of 15th International Symposium on Remote 
Sensing of Environment .Ann Arbor, MI, 125-144. 

Barringer, T., V.B. Robinson, J.C. Coiner, and R.C. Bruce. 1980. 
"Landsat Analysis of Tropical Rain Forest Succession Employing a 
Terrain Model", Proceedings of 14th International Symposium on Remote 
Sensing of Environment . San Jose, Costa Rica, 1691-1700. 

Bartels, C.P.A. and R.H. Ketellapper, 1979, Exploratory 
Statistical Analysis of Spatial Data , Martinus Nijhoff , Boston. 

Bartlett, M.S., 1969, "Two-Dimensional Nearest-Neighbor Systems 
and Their Ecological Application", Proc. Inti. Symp. Ecol. , New Haven, 
Conn. 


, 1975, The Statistical Analysis of Spatial Pattern . 

Chapman and Hall, London. 

Bassett, K., 1972, "Numerical Methods for Map Analysis", Prog. 

Geog. . Vol 4. 

Batty, M. and R. Sammons, 1978, "On Searching for the Most 
Informative Spatial Pattern" , Envt. and Planning A. 

Beeftink, W.G., 1979, "Vegetation Dynamics in Retrospect and 
Prospect Introduction to the Proceedings of the Second Symposium of the 
Working Group on Succession Research on Permanent Plots", 
Vegetatio . 40(2): 101 - 105. 

Bell, E.J., 1974, "Markov Analysis of Land Use Change: An 

Application of Stochastic Processes to Remote Sensing Data", Socio- 
Economic Ping. Sci. , (8):311 - 316. 

,and R. Hinojosa. 1977, "Markov Analysis of Land Use 

Change: Continuous Time and Stationary Processes", Socio-Economic Ping. 
Sci., 11:13-17. 


42 


tollhouse, D.R. 1977, "Some Optimal Designs for Sampling in Two 
Dimensions", Biometrika , 64(3): 60S - 611. 

Bennett, R.J. 1979, Spatial Time Series . Pion, London. 

and R.J. Chorley, Environmental Systems: Analysis, 

Forecasting and Control , Princeton . 

and R.J. Chorley, 1978, Environmental Systems: 

Philosophy, Analysis and Control , Methujn, London. 

Berry, B.J.L., "Sampling Coding and Storing Flood Plain Data", 
USDA, Farm Econ. Div., Agric. Handbk. 137. 

tosag, J. 1974, "Spatial Interaf tion and the Spatial Analysis of 
Lattice Systems", J. Roy. Stat. Sor., 3, 36:192 - 235. 

and P.J. Digglf 1977, "Simple Monte Carlo Tests for 

Spatial Pattern" , Appl. Stat. . ?5(3):327 - 333. 

Box, G.E.P. and G.M. Je.iKins, 1976, Time Series Analysis: 
Forecasting and Control . Rev. Ed., Holden - Day, San Francisco. 

Brandsma, A.S. and R.H. Ketellapper. 1979, "Further Evidence on 
Alternative Proceedures for Testing of Spatial Autocorrelation Among 
Regression Disturbances" in Bartels and Ketellapper. 1979. 

Campbell, J.B. 1981, "Spatial Correlation Effects Upon Accuracy 
of Supervised Classification of Land Cover", Phot. ENG. Rem. Sens. . 
47(3) : 355 - 363. 

Clements, F.E., 1916, "Plant Succession: An Analysis of the 

Development of Vegetation" , Carnegie Inst. Pub. No. 242 , Washington . 

. 1936, "Nature and Structure of the Climax", Jrnl. 

Ecology . 24:552 - 584. 

Cliff, A.D. 1977, "Quantitative Methods: Time Series Methods for 
Modelling and Forecasting", Prog. Hum. Geog. , (1):492 - 502. 

. 1980, "Quantitative Methods: Spatial Modelling", 

Prog. Hum. GeoR. , (4): 568 - 576. 

, P. Haggett, J.K. Ord, K. Bassett, R. Davies. 1975. 

Elements of Spatial Structure , Cambridge Univ. , Cambridge. 

. and J.K. Ord. 1975a, "The Comparison of Means When 
Samples Consist of Spatially Autocorrelated Observations", Envt. & 
Ping.. (7): 725 - 734 


43 



. 1975b, "Model Building and the Analysis 

of Spatial Pattern in Human Geography",^ Roy. Stat. Soc. , ji, (37): 
297-348. 

Coiner, J.C., 1980, "Using Land sat to Monitor Changes in 

Vegetation Cover Induced by Desertification Processes", Proc. 14th 
Inti. Symp. Rea. Sens. Envt. , Envtl. Res. Inst, of Mich., Ann Arbor. 

Connell, J.H. and R.O. Slatyer. 1977, "Mechanisms of Succession 
in Natural Communities and Their Role in Community Stability and 
Organization", Am. Nat. . 111(982): 1119 - 1143. 

Craig, R. G. 1979, "Autocorrelation in Landsat Data”, Proc. 13th 
Inti. Symp. Rem. Sens. Envt. , Envtl. Res. Inst. Mich., Ann Arbor. 

, and M. L. Labovitz. 1980, "Sources of Variation in 

Landsat Autocorrelation" , Proc. 14th Inti. Symp. Rem. Sens. Envt. , 
Envt. Inst. Mich., Ann Arbor. 

Curtis, J.T., and R.P. McIntosh. 1950, "The Interrelationships of 
Certain Analytic and Synthetic Phytosociological Characters", Ecology , 
31(3):434 - 455. 

Drury, W.H. and I.C.T. Nisbet. 1973, "Succession", J. Arnold 
Arb. . 54(3):331 - 368. 

Dutta, M. 1975, Econometric Methods . Southwestern Publishing 
Co. .Cincinnati. 

Estes, J.E., J.R. Jensen and D.S. Simonett. 1980, Rem. Sens. 
Envt. , 10(1 ):43 - 80. ' 

Fitzpatrick-Lins, K. 1981, "Comprison of Sampling Proceedures and 
Data Analysis for a Land-Use and Land -Cover Map", Phot. Eng. Rem. 
Sens. . 47(3):343 - 351. 

Getls, A. and B.N. Boots. 1978, Models of Spatial Processes . 
Cambridge Univ., Cambridge. 

Gleason, H.A. 1926, "The Individualistic Concept of Plant 
Association", Torrey Bot. Club Bull. , (53): 7 - 26. 

Gould, P.R. 1970, "Is Statistix Inferens the Geographical name 
for a Wild Goose?", Econ . Geog . , 46 (Supp.):439 - 484. 

Granger, C. W. 1969, "Spatial Data and Time Series", Stud. Reg. 
Sci. , A. J. Scott (Ed.), Pion, London. 


44 


Greig-Smith, P. 1964, Quantitative Plant Ecology , 2nd ed., 
Butterwor th , London . 

Haggett, P. 1963, "Regional and Local Components in Land -Use 
Sampling", Erdkunde , 7 : 108 - 114. 

Haggett, P., A.D. Cliff and A. Frey. 1977, Locational Models in 
Human Geography . 2nd ed., Halsted Press, New York. 

Haining, R. P. 1978, "Specification and Estimation Problems in 
Models of Spatial Dependence", Studies in Geog. #24 , Dept. ofGeog., 
Northwestern Univ., Evanston. 

Hilborn, R. 1979, "Some Long Term Dynamics of Predator-Prey 
Models With Diffusion", Ecol. Mod. , (6):23 - 30. 

van Hulst, R.V. 1978, "On the Dynamics of Vegetation: Patterns of 
Environmental and Vegetational Change", Vegetatio , 38(2): 65 - 75. 

. 1979a, "On the Dynamics of Vegetation: Succession in 

Model Communities", Vegetatio , 39(2):85 - 96. 

. 1979b, "On the Dynamics of Vegetation: Markov Chains 
as Models of Succession", Vegetatio 40(1): 3 - 14. 

Jensen, J.R. 1979, "Spectral and Textural Features to Classify 
Elusive Land Cover at the Urban Fringe", Prof. Geog. , 31(4):400 - 408. 

Kooijman, S.A.L.M. 1980, "Inference About Dispersal Patterns", 
Acta BioTheor , 28(3): 149 - 189. 

Krumbein, W.C. and F.A. Graybill. 1965 . Introduction to 
Statistical Models in Geology , McGraw-Hill : New York . 

Lebart, L. 1969, "Analyse Statistique de la Contiguite", Publ. 
Inst. Stat. Univ. Paris , 18:81 - 110. 

Levin, S.A. and R.T. Paine, "Disturbance, Patch Formation, and 
Community Structure", Proc. Nat. Acad. Sci., USA , 71(7): 2744-2747. 

Lillesand, T. M. and R. W. Kiefer. 1979, Remote Sensing and Image 
Interpretation , Wiley, New York. 

Martin R.L., and J.E. Oeppen. 1975, "The Identification of 
Regional Forecasting Models Using Space-Time Correlation Functions", 
Inst. Brit. Geog. , 66:95 - 118. 

Mather, P.M. 1976, Computational Methods of Multivariate Analysis 
in Physical Geography , John Wiley & Sons : New York. 


45 



T v T 







Miles, J. 1979, Vegetation Dynamics , Chapman and Hall, London 

Miller, L.D., K. Nualchawee and C. Tom. 1979, "Analysis of the 
Dynamics of Shifting Cultivation in the Tropical Forests of Northern 
Tnailand Using Landscape Modeling and Classification of LANDSAT 
Imagery". Proc. 12th Inti. Symp. Rem. Sens. Envt^, Envtl. Res. Inst. 
Mich. 

Odum, E.P. 1969, "The Strategy of Ecosystem Development", 
Science, 164:262 - 270 

Omakupt, M. and V. Vunpiyarat. 1980, "Remote Sensing for Land Use 
Inventory", Proc. 14th Inti. Symp. Rem. Sens. Envt., Envtl. , Res. Inst. 
Mich., Ann Arbor. 

Pielou, E.C. 1969, An Introduction to Mathematical Ecology . 
Wiley, New York. 

Reeves, R.G., (Ed.). 1975, Manual of Remote Sensing( 2 vols.). Am. 
Soc. Photogram., Falls Church. 

Robinson, V.B. 1980, "On the Use of Markovian Equilibrium Vectors 
in the Evaluation of Land Use Policy" , Soc io-Economlc Planning Sciences , 
14, 85-89. 

Robinson, V.B. 1978. "Information Theory and Sequences of Land 
Use: An Application", The Professional Geographer , 30 : 174-179. 

Robinson, V.B., T. Barringer, J.C. Coiner. 1982. "Dynamic 
Modeling of Vegetation Change in Arid Lands", Proceedings of First 
International Thematic Conference on Remote Sensing , 121-131 , Cairo, 

Egypu 


Rogers, A. 1974, Statistical Analysi s of Spatial Dispersion . 
Pion, London. 


1 


( 


Shugart, H.H., Jr., M.S. Hopkins, I.P. Burgess, and A.T. 
Mortlock, 1980, "The Development of a Succession Model for Subtropical 
Rain Forest and its Application to Assess the Effects of Timber Harvest 
at Wiangaree State Forest, New South W*ies", JL Envt . Mgt . , 11:243 - 
265. 

Siatyer, R.O. (Ed.). 1977, "Dynamic Changes in Terrestrial 
Ecosystems: Patterns of Change, Techniques for Study and Applications 
to Management", UNESCO . 

Smith, J.M. 1974, Models in Ecology , Cambridge Ur.iv , Cambridge. 

Smith, T.E. 1980, "A Central Limit Theorem for Spatial 
Samples". Geog. Anal. , 12(4):299 - 32*. 


46 



Stow, D.A., L.R. Tinnev, and J.E. Estes, "Deriving Land Use/Land 
Cover Change Statistics from Landsat: A Study of Prime Agricultural 

Land", Proc . 14th Inti. Svmp. Rem. Sena. Envt. . Envtl. Rea. Inst. 
Mich., Ann Arbor. 

Strahler, A.H. and C.E. Woodcock, "Stratification of Forest 
Vegetation for Timber Inventory Using Landsat Digital Imagery and 
Collateral Data", Proc. 14th Inti. Symp. Rem. Sens. Envt. . Envtl. Res. 
Inst. Mich., Ann Arbor. 

Swain, P.H. and S.M. Davis. 1978, Remote Sensing: The 

Quantitative Approach , McGraw-Hill, New York. 

Taranik, J.V. 1978, "Principles of Computer Processing of Landsat 
Data for Geologic Applications " , U.S.G.S. , Open File Report 78-117 , 
Sioux Falls, South Dakota. 

Tobler, W.R. 1967, "Of Maps and Matrices", J. Reg. Sci. , 7(2):275 

- 280. 


. 1970, "A Computer Movie Simulating Urban Growth in 

the Detroit Region", Econ . Geog . , 46:234 - 240. 

Tubbs, J.D. and W.A. Coberly,1978, "Special (sic) Correlation and 
Its Effect Upon Classification Results in Landsat", Proc. 12th Inti. 
Symp. Rem. Sens. Envt. , Envtl . Res . Inst . Mich . , Ann Arbor. 

Whittaker, R.H. 1953, "A Consideration of Climax Theory: The 

Climax as a Population and Pattern", Ecol . Monogr . , 23:41 - 78. 

Whittle, P. 1954, "On Stationary Processes in the Plane", 
Biometrika , 41:434 - 449. 

. 1963, "Stochastic Processes in Several Dimensions", 

Bull. Inti. Stat. Inst ,. 34:974 - 993. 

Woodcock, C.E., A.H. Strahler and T.L. Logan. 1980, 

"Stratification of Forest Vegetation for Timber Inventory Using Landsat 
and Collateral Data", Proc. 14th Inti. Symp. Rem. Sens. Envt. . Envtl. 
Res. Inst. Mich., Ann Arbor 

Yarranton, G.A. and R.G. Morrison. 1974, "Spatial Dynamics of a 
Primary Succession: Nucleation", Jrnl. Ecol. , 62:417 - 428. 

Yule, G.U. and M.G. Kendall. 1957, An Introduction to the Theory 
of Statistics, Griffin, London. 


47 


